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Abstract 

Our  ab  initio  theoretical  calculations  indicate  that  the  LiFLi' 
anion  possesses  two  bound  states  which  are  of  and 

symmetry  and  which  have  vertical  detachment  energies  of  0.94  and 
0.89  eV,  respectively.  The  LiFLi+  cation  and  the  first  four  electronic 
states  of  the  neutral  LiFLi  were  also  considered.  The  Renner-Teller 
splitting  for  the  excited  state  of  the  neutral  is  reported.  The 

neutral  and  both  ions  are  thermodynamically  stable. 
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I.  Introduction 


Atomic  and  molecular  anions  usually  do  not  support  more  than 
one  bound  electronic  state,  but  exceptions  have  been  found  in 
electron  photodetachment  experiments.  Atoms  with  two  or  more 
bound  anionic  electronic  states  are  covered  in  Ref.  [1]  and  molecular 
cases  are  listed  in  Ref.  [2]. 

Here  we  report  results  of  our  ab  initio  calculations  on  two 
electronically  and  geometrically  stable  states  of  LiFLi'.  This  species 
may  be  considered  as  an  example  of  the  lithium  substituted  double- 
Rydberg  (DR)  anion  HFH*  studied  earlier  [3-7].  DR  anions,  such  as  the 
Td  isomer  of  NH4'  or  the  Csv  isomer  of  HaO’,  consist  of  a  closed-shell 
cation  core  plus  two  extra  electrons  described  by  orbital(s)  of 
Rydberg  nature  [6,7], 

Ha*,  which  is  an  obvious  candidate  for  a  DR  anion,  is 
electronically  unstable  at  the  Dah  equilibrium  geometry  of  Ha"^  [3,4]. 
In  contrast,  we  recently  found  that  the  analog  alkali  metal  trimers 
Lia*  and  Naa*  support  more  than  one  geometrically  and  electronically 
bound  anionic  state  [8].  The  linear  ground  l2g+  states  of  these  alkali 
metal  trimer  anions  were  probed  in  electron  photodetachment 
experiments  [9],  but  triplet  and  quintet  bound  states  of  these  species 
characterized  in  Ref.  [8]  have  not  been  observed  yet. 

Encouraging  results  on  the  alkali  metal  trimers  prompted  us  to 
study  the  lithium  substituted  DR  anions  and  the  underlying  neutral 
and  cation  species.  In  this  work,  we  report  on  LiFLi  and  its  cation 
which  have  recently  been  theoretically  studied  by  Rehm  et  al.  [10]  in 
their  extensive  work  on  superalkalis  (i.e.,  molecular  systems  with 

first  ionization  potential  smaller  than  that  of  the  Cs  atom).  To  the 

best  of  our  knowledge,  the  LiFLi*  anion  has  not  been  experimentally 
studied  yet. 

II.  Computational  Aspects 

For  the  lithium  atom,  we  used  the  Dunning  (9s5p/3s2p)  one- 

electron  basis  set  [11]  supplemented  with  diffuse  s  and  p  functions 
with  the  same  exponent  0.0074  [12]  and  one  d  function  with  the 


exponent  0.2  [13].  This  basis  set  is  detailed  in  Ref.  [8].  For  the  oxygen 
atom  we  employed  Dunning's  aug-cc-pVDZ  basis  set  which  was 
designed  to  describe  anionic  species  [14J.  Cartesian  d  functions  were 
used  throughout  the  calculations  and  the  full  basis  set  for  the 
triatomic  consists  of  63  contracted  gaussian  functions. 

Potential  energy  surfaces  were  explored  within  the  complete 
active  space  (CAS)  self-consistent  field  (SCF)  formalism.  Molecular 
orbitals  which  result  from  the  core  Is  atomic  orbitals  were  doubly 
occupied  in  every  configuration  state  function  (CSF).  The  neglected 
core-core  and  core-valence  correlation  effects  are  negligible  for  the 
lithium  and  fluorine  atoms  due  to  the  low  polarizability  of  the  Is 
cores.  The  remaining  eight  (cation),  nine  (neutral)  or  ten  (anion) 
electrons  were  distributed  in  all  possible  ways  among  four  ai,  two  bi, 
and  two  b2  molecular  orbitals  (the  C2v  symmetry  labeling  is  used). 
This  choice  of  the  active  space  led  to  1764,  2352,  1176,  and  1512 
CSF’s  for  the  cation,  neutral,  singlet  and  triplet  anion,  respectively. 
The  CAS  SCF  calculations  were  performed  with  the  Utah  MESS-KIT 
modular  electronic  structure  codes  [15]  which  generate  analytical 
second  geometrical  derivatives.  Stationary  points  on  the  potential 
energy  surfaces  were  determined  using  our  automated  surface 
walking  algorithm  [16]. 

The  restricted  CAS  SCF  approach  is  capable  to  predict  accurate 
geometries  but  it  is  inappropriate  to  accurately  compare  energies  of 
species  with  a  different  number  of  electrons.  Hence  we  employed  the 
quadratic  configuration  interaction  (QCI)  approach  with  single, 
double  and  approximate  triple  excitations  (SD(T))  [17]  to  determine 
the  final  detachment  energies  and  electron  affinities.  The  QCISD(T) 
approach  is  size-extensive  and  takes  into  account  dynamical 
correlation  affects.  In  the  QCI  calculations,  we  kept  the  six  core 
electrons  uncorrelated.  We  checked  that  this  restriction  changes  the 
vertical  detachment  energy  for  the  ground  state  anion  by  less  than 
0.003  eV.  The  QCI  results  were  obtained  with  the  Gaussian  92  suite 
of  codes  [18]. 


III.  Results 
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The  stationary  points  on  the  potential  energy  surfaces  of  the 
cation,  neutral,  and  anion  are  characterized  at  Table  I.  me  harmonic 
vibrational  frequencies  were  determined  at  the  CAS  SCF  level  using 
analytical  second  derivative  methods.  The  spatial  extents  of  the 
electronic  charge  distributions  are  characterized  by  the  SCF  values  of 
<R2>  in  column  four.  The  values  of  the  relative  energies  reported  in 
the  fifth  column  were  obtained  at  the  QCISD(T)  level.  The  energies  of 
the  first  four  electronic  states  as  functions  of  the  bending  angle  are 
displayed  in  Fig.  1.  These  energies  were  obtained  in  a  state  averaged 
(SA)  CAS  SCF  calculations  [19]  in  which  the  orbitals  were  optimized 
simultaneously  for  the  four  electronic  states.  Our  vertical  (VDE  & 
VIP)  and  adiabatic  (ADE  &  AIP)  detachment  and  ionization  energies 
are  presented  in  Table  II. 

A.  The  Cation  and  Neutral 

Our  geometries  for  the  closed-shell  cation  and  the  ground  state 
neutral  are  in  a  good  agreement  with  the  results  of  Rehm  et  al.  [10]. 
The  cation  has  a  linear  structure,  whereas  the  ground  state  of  the 
neutral  is  bent  (^Ai)  with  a  barrier  for  linearization  of  5.3  kcal/mol. 
The  neutral's  5ai  highest  occupied  molecular  orbital  (HOMO)  is 
dominated  by  the  Li's  2s  orbitals  which  interact  constructively  with 
each  other  and  destructively  with  small  s-  and  Pz-type  contributions 
from  the  central  F  atom.  The  constructive  interaction  between  the  Li 
ligands  is  apparently  able  to  overwhelm  the  repulsion  between  the 
positive  charges  that  accumulate  on  these  centers  and  cause  the 
species  to  be  bent.  In  the  neutral's  first  excited  state,  however,  the 
unpaired  electron  occupies  the  3b2  orbital  in  which  the  Li's  2s 
orbitals  interact  destructively.  Hence,  the  equilibrium  structure  of 
the  first  excited  state  of  LiFLi  is  linear  (2Zu+). 

For  the  (^Hu)  next  excited  state  of  the  neutral,  the  equilibrium 
structure  is  again  linear.  The  corresponding  2ku  HOMO  is  dominated 
by  the  2p  Li  orbitals.  Contributions  from  the  F’s  p-type  orbitals  are 
much  smaller  and  the  Li-F  p-p  interaction  is  destructive.  Due  to 
second-order  Renner-Teller  coupling,  the  bending  frequencies  are 
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not  degenerate,  and  distortion  along  the  lower  frequency  mode  leads 
to  the  2b  1  state  in  C2v  geometry. 

Our  AIP's  for  the  2Ai,2iy+,  and  2nu  states  of  LiFLi  are  3.93, 
3.31,  and  2.29  eV,  respectively.  The  experimental  ionization  energy 
for  LiFLi  is  not  available,  but  the  result  for  NaClNa  of  4.15  ±  0.2  eV 
[20]  is  quite  similar  to  our  finding. 

Let  us  recall  that  the  cation  HFH+  has  a  C2v  structure  (the  angle 
is  1130),  and  that  the  ground  2a  i  state  of  the  HFH  radical  is 
geometrically  unstable  with  respect  the  b2  vibrational  mode.  The 
excited  2b 2  state  of  HFH  has  a  C2v  geometry  with  bond  lengths  quite 
different  from  the  cation’s  equilibrium  structure.  Finally,  the 
corresponding  2b  |  state  of  HFH  has  a  geometry  close  to  the  cation 
equilibrium  C2v  geometry  [4]  Clearly,  the  substitution  of  H’s  by  Li's 
leads  to  remarkable  changes  in  the  equilibrium  geometries  of  the 
cation  and  neutral  electronic  states.  Most  of  these  differences  can  be 
attributed  to  (i)  the  fact  that  Li  is  much  more  electropositive  than  H, 
and  (ii)  Li  has  2s  and  2p  valence  orbitals  to  use  but  H  has  only  Is 
valence  orbitals. 

B.  The  Anion 

The  ground  electronic  state  of  LiFLi*  has  a  linear  equilibrium 
structure  with  a  geometry  close  to  that  of  LiFLi+.  The  wavefunction 
of  the  anionic  state  can  be  viewed  as  that  of  LiFLi+  with  a 

single-  coupled  electron  pair  "attached",  one  electron  residing  on  each 
of  the  two  Li  centers.  This  wavefunction  appears  in  our  calcualtions 
as  a  two-configurational  function  with  the  Cl  mixing  coefficients  of 
the  ag2  and  Ou2  configurations  equal  to  0.76  and  -0.62,  respectively. 
The  Og  and  Ou  orbitals  are  dominated  by  Li's  s-type  basis  functions. 
The  bonding  interaction  between  ligands  resulting  from  the  Og2 
contribution  is  counterbalanced  by  the  antibonding  interaction 
coming  from  the  Ou2  contribution  In  addition,  the  Coulomb  repulsion 
between  effective  negative  charges  on  every  atom  favors  the  linear 
structure.  The  Og  and  Ou  orbitals  are  much  more  diffuse  than  in  the 
ground  ^Ig"^  or  first  excited  2i„+  state  of  the  neutral  radical  as 
reflected  by  the  values  of  <R2>  in  Table  I. 
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For  the  ground  state  of  the  anion,  the  VDE  to  the  state  of 

the  radical  is  0.94  eV,  which  aufers  by  0.23  eV  from  the  ADE,  (i.e.,  by 
the  energy  of  linearization  of  the  i  radical).  The  anion  is  also 
electronically  stable  at  the  C2v  equilibrium  geometry  of  the  i 
radical,  but  only  by  0.60  eV. 

In  addition  to  the  ground  state,  the  anion  LiFLi-  possesses 

another  electronically  bound  state,  again  with  an  equilibrium 

geometry  close  to  that  of  the  cation  and  the  ground  state  anion.  This 
function  can  be  viewed  as  a  LiFLi"*"  cation  with  two  "attached" 
electrons,  one  in  a  Og  orbital  and  the  second  in  the  corresponding  Ou 
orbital.  The  Og  orbital  is  as  described  above  and  the  Ou  orbital  is 
dominated  by  the  Li's  s-  and  pz-type  basis  functions,  whereas  the 
contribution  from  the  F's  pz  functions  is  much  smaller.  The  s-pzLiF 
interactions  are  antibonding,  and  the  pz-pz  are  bonding.  The  SCF 
values  of  <R2>  predicts  to  b®  loss  extended  than  the  ground  ^Zg"*" 
state,  which  may,  however,  be  an  artifact  of  the  one-configuration 
approximation  used  in  calculating  <R2>.  The  VDE's  from  the  anion's 
^Zu"*"  state  to  the  2Zg+  and  2Zu+  states  of  the  radical  are  0.89  and  1.28 
eV,  respectively.  The  triplet  anion  is  also  electronically  stable  at  the 
2Ai  minimum  of  the  radical,  but  only  by  0.38  eV.  Hence,  we  conclude 
that  the  anion  ^Zu"*"  is  electronically  and  geometrically  stable. 

All  of  the  species  considered  here  are  thermodynamically 
stable  at  the  QCISD(T)  level.  The  cation  is  thermodynamically  stable 
with  respect  to  Li'*‘0S)  +  LiF  by  63.0  kcal/mol  (with  zero-point 
vibrational  corrections  included).  The  neutral  in  the  ground  state  is 
stable  with  respect  to  Li(2S)  +  LiF  by  30.8  kcal/mol,  and  the  anion  in 
its  ground  state  is  stable  by  33,5  kcal/mol  with  respect  to  Li*(^S)  + 
LiF.  The  larger  stability  of  the  anion  compared  to  that  of  the  neutral 
reflects  the  fact  that  electron  affinity  is  larger  for  LiFLi  than  for  Li. 

IV.  Conclusions 

Theoretical  calculations  indicate  that  lithium  substituted 
double-Rydberg  anions  can  possess  two  bound  anionic  states.  The 
fully  symmetric  singlet  state  is  the  ground  state.  In  addition,  the 
triplet  state  is  electronically  bound.  The  analog  triplet  states  of  Li2F-, 
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Li30-,  LUN"  have  symmetry  35:u+,  3e',  and  3X2,  respectively.  Our 
results  for  Li2F’  were  reported  here;  for  LisO*  and  LuN",  the  triplet 
state  is  subject  to  first-order  Jahn-Teller  distortion  and  complete 
results  will  be  discussed  elsewhere. 

In  the  case  of  Li2F',  the  vertical  and  adiabatic  detachment 
energies  for  the  ground  state  anion  are  0.94  and  0.72  eV, 

respectively.  The  ground  state  of  the  radical  has  a  Czv  equilibrium 
geometry  with  the  linearization  barrier  (0.23  eV)  equal  roughly  to 
the  difference  between  the  two  above  energies. 

For  the  3i„+  anion  the  vertical  detachment  energies  to  the 
and  25;^+  neutral  states  are  0.89  and  1.28  eV,  respectively.  The 
corresponding  adiabatic  energies  are  0.66  and  1.28  eV.  Both  the 
singlet  and  triplet  anionic  state  remains  electronically  stable  in  the 
C2v  equilibrium  geometry  of  the  ground  state  radical. 

Excited  electronic  states  of  LiFLi  have  never  been  studied 
before  experimentally.  Both  ^1^+  and  ^riu  states  have  minima  close  to 
the  cation  equilibrium  geometry.  They  are  separated  from  the 
ground  state  by  0,63  and  1.64  eV,  respectively.  The  bending 

frequencies  in  the  ^riu  state  split  by  23  cm*^  due  to  the  second-order 
Renner-Teller  coupling. 

The  cationic,  neutral,  and  anionic  LiFLi  species  should  be 
amenable  to  experimental  studies  because  they  are 
thermodynamically  stable  to  dissociation. 
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Table  I.  Geometries  (distances  in  bohr.angles  in  degrees)  and  CAS  SCF 
harmonic  frequencies  (cm-1)  for  the  cation,  neutral,  and  anion  of 
LiFLi-  Spatial  extent  of  the  SCF  electronic  charge  distribution  <R2>  is  in 
(bohr)2,  and  the  relative  QCISD(T)  energies  are  in  kcal/mol. 


Species/State 

Geometry 

Frequencies 

<3l2> 

Energy* 

LiFLi+  lSg+ 

Dooh 

RLiF=3.209 

151(jtu);614(Og); 

767(Ou) 

58 

0.0 

LiFLi  2a  1 

C2v 

RLiF=3.235 

Z  =97.94 

249(ai);643(b2); 

675(ai) 

86 

-90.7 

LiFLi  2£g+ 

Dooh 

RLiF=3.172 

143i(jCu);627(og); 

565(au) 

103 

-85.4 

LiFLi  2Xy+ 

Dooh 

RLiF=3.181 

119(jCu);627(og); 

894(Ou) 

1  10 

-76.2 

LiFLi  2nu 

Dooh 

RLiF=3.157 

71(7Cu->^Bi); 

94(71u->2Ai); 

643(ag);  751(Ou) 

279 

-52.9 

LiFLi-  lZg+ 

Dooh 

RLiF=3.197 

132(3Cu);598(ag); 

752(<Tu) 

224 

-107.2 

LiFLi-  3iu+ 

E>ooh 

RLiF=3.203 

135(jcu);595(ag); 

745(au) 

187 

-105.9 

^Electronic  energies  without  zero-point  vibrational  correction 
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Table  11.  Vertical  and  adiabatic  detachment  energies  (VDE  &  ADE) 
and  adiabatic  ionization  potentials  (AIP)  (in  eV)  for  the  anionic  and 
neutral  LiFLi. 


Species 

Transition 

Detachment  or 
Ionization  Energy 

LiFLi- 

lZg+  ^  2i:g+;  2Ai 

0.94  (VDE);  0.72  (ADE) 

LiFLi- 

3Xu+  ->  2ig+;  2a  1 

0.89  (VDE);  0.66  (ADE) 

LiFLi* 

3Xu+ 

1.28  (VDE);  1.28  (ADE) 

LiFLi 

2Ai  ^  llg+ 

3.93  (AIP) 

LiFLi 

2Xu+  ^  llg+ 

3.30  (AIP) 

LiFLi 

2nu+  ^  1 V 

2.29  (AIP) 

Figure  Caption 
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Figure  1.  SA  CAS  SCF  energies  of  the  first  four  electronic  states  of 
LiFLi  with  respect  to  the  bending  angle.  An  angle  of  180  deg 
corresponds  to  the  linear  structure. 


Energy  (kcal/mol) 


